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Abstract 

We present a second-quantized field theory of massive spin one-half particles or antiparticles 
in the presence of a weak gravitational field treated as a spin two external field in a flat Minkowski 
background. We solve the difficulties which arise from the derivative coupling and we are able 
to introduce an interaction picture. We derive expressions for the scattering amplitude and for 
the outgoing spinor to first-order. In several appendices, the link with the canonical approach in 
General Relativity is established and a generalized stationary phase method is used to calculate 
the outgoing spinor. We show how our expressions can be used to calculate and discuss phase 
shifts in the context of matter-wave interferometry (especially atom or antiatom interferometry). 
In this way, many effects are introduced in a unified relativistic framework, including spin- 
gravitation terms: gravitational red shift, Thomas precession, Sagnac effect, spin-rotation effect, 
orbital and spin Lense-Thirring effects, de Sitter geodetic precession and finally the effect of 
gravitational waves. A new analogy with the electromagnetic interaction is pointed out. 



1 Introduction 

The development of high accuracy atom interferometers, used as clocks in the microwave or in the 
optical domain, as inertial sensors (gyros, gravimeters, gradiometers...) or for the determination 
of atomic masses and of the fine structure constant ^, ^, ^, ^, ^, |^, |^, |9|, requires 
now a framework to describe the interference of atom waves in a rigorous way. On one hand, 
one needs to investigate general relativistic effects including those involving the spin of the atoms 
and, on the other hand, it is necessary to take into account the statistical properties (bosonic or 
fermionic) of the interfering particles, given the development of coherent atom wave sources and 
also for a proper treatment of the detection noise. One would also like to be able to discuss the 
propagation of antimatter in interferometers in the presence of gravitation and as suggested in 



reference [22| the properties of coherent antimatter waves (generated by an antiatom laser such 
as an antihydrogen Bose-Einstein condensate). This is possible only within relativistic quantum 
field theory. Atoms (or antiatoms), in a given eigenstate of the internal atomic Hamiltonian, are 
considered as elementary particles having a rest mass fixed by the energy of the atomic level and a 
spin equal to the total angular momentum of the atom in that level. In this paper, we shall consider 
only Dirac particles for illustration. This is the simplest example of particles with spin which still 
contains most interesting effects related to spin and applies to neutron or electron interferometry as 
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a special case. The generalization to others spin values is possible along similar lines with Dirac-type 
equations (Bargmann-Wigner [^], de Broglie fusion method |28|, Durand [^^). Our overall goal 



is to introduce gravitation and general relativistic effects at the quantum level of modern atomic 
physics and quantum optics experiments. For this, we propose an extension of our first paper on 
atom interferometry in General Relativity which includes now a second-quantization scheme 
for the atom waves in the presence of gravitational and electromagnetic fields. The point of view 
adopted in this paper is the extrinsic point of view using purely quantum field theory in a flat 
Minkowski background. The connection with the canonical intrinsic approach, using Dirac equation 
in curved space-time, is made in Appendices A and B. The reader who wishes to start with this 
canonical approach is thus invited to read first these appendices. 

In the main text, we begin right away with the minimal coupling Lagrangian in flat space-time, 
which is identical to the one derived in curved space-time for a standard choice of tetrads. Then, 
we proceed with the quantization of the Dirac fleld and we emphasize the difficulties which arise 
because of the derivative coupling. These difficulties are solved in a consistent scheme which allows 
also to define an interaction picture. The evolution operator and the S'-matrix are constructed 
and we demonstrate explicitly a conjecture of Gupta. These results are used to derive formulas 
for the transition amplitude and for the outgoing spinor in the weak-field approximation, first in 
configuration space and second in the momentum representation. An expansion in the perturbation 
wave vector (hk/mc) is used to retrieve various physical effects, some of which are well-known. A 
new analogy with the electromagnetic interaction is presented which includes all components of the 
field h^'^ and generalizes gravitoelectric and gravitomagnetic interactions. In Appendices C and D, 
calculations of the outgoing spinor are sketched, first with a generalized stationary phase method in 
configuration space and second in the momentum representation. 

2 Lagrangian theory 

Considered as a field theory in fiat spacetime, the theory describing the interaction of matter with 
a given gravitational field will be defined by a Lagrangian density of the following form |15, 19|: 



C = Co- Ih^'^T^, , (1) 



where h^'^ is the given external field, and where Cq is the free Lagrangian density of the matter 
field and T^i, the corresponding stress-energy tensor. For a Dirac field, one has respectively]^, in 
symmetrical form. 



he 



(2) 
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The total Lagrangian density (|l|) then becomes 

£ = (1 + i/i) £o - lihchf^" ^ (tm^- - "^-7^) (4) 



where^ h = h'^^ = rj^^yh^'^ . This Lagrangian density can also be considered as obtained from the 
Lagrangian density valid in General Relativity for the interaction of the Dirac field with a prescribed 
gravitational field in the linear approximation (see the Appendices). 

'^The conventions used here for the metric and for the Dirac equation and matrices are generaUy those of [p^. In 
particular, the signature of the metric is taken as (+, — , — , — ). Greek indices /i, i^, ... run from to 3 and latin indices 
run from 1 to 3. The space-time 4-vector is written x = (a;", x) = {ct, x). The partial derivatives with the right arrow 
act on the right and those with the left arrow on the left. 

^The term {h/2)£o in comes from the term containing Co in and was omitted in jl6j. However, Co which 
vanishes when the free Dirac equation is satisfied, does not vanish here. 



The equations for ^' and derived from (Q) are the following: 



[(1 + i/i) (^7^^ - mc/h) - '-h^^^^^^t, - ^a./i'^^T^ + ^9^/^71^ = 0, (5) 

^[(-^7'^'^^ - mc/h) (1 + Ih) + '-"d.j^h^^'' + '-d^h'^'^j^ - '-d^h-f^'] = 0. (6) 

They admit the conserved current 

f = c^[^/^ + 1/17/^ - ^h^''^,]^, (7) 

which coincides with the usual Dirac current when the gravitational field vanishes. 

In order to stress some other differences and analogies with electromagnetism, we may write the 
equation for ^ with a covariant derivative in the usual sense of non-Abelian gauge field theories in 
flat space-time 



where 



<y^'' = \{i^Y-Yin (9) 

and where the factor (l + i/i) has been removed. The Poincare generators are associated with gauge 
fields and a local gauge invarianceQ 



1 + {dxi,-d,ix)-edx 



■^{x) (10) 



h'^j^A^) = h^y{x) - df^iy - dyif, . (11) 

3 Quantization 

We want now to proceed with the quantization of the field \I' submitted to the interaction defined 
by (§). When they are applied directly, the standard methods lead to some difficulties coming from 
the presence of a derivative coupling in the Lagrangian. These problems are, first, briefly discussed, 
then, a solution is presented allowing the quantization together with the definition of the interaction 
picture. 



3.1 Difficulties with the derivative coupling 

Two methods can be used, a priori, according to whether the Lagrangian is taken under a symmet- 
rical or a asymmetrical form. 



3.1.1 Symmetrical Lagrangian. 

Starting from the symmetrical Lagrangian (^), the usual anticommutation relations will be obtained 
from the following expression of the conjugate momentum 

dC 

n = 2^, with ^! = dt^, (12) 

^If one considers these equations as first-order equations with respect to the h'^^'s, the factor (1 + h/2) in the first 
term can be replaced by 1. But it can be shown that the equations so obtained cannot be derived from a Lagrangian, 
if this latter is restricted to depend iinearly on ^ and and to admit first-order derivatives oniy. 

^Differences with usuai Yang-Miils theories come from the non-commutation of the Lorentz generators with the 
Dirac matrices and from the fact that the translation generators act on space-time itself. 



which gives 
and 



(13) 
(14) 



{^'a(x), n^(y)}^o=j^o = ihSaiB S{x - y) 

where a and /? are spinorial indices. 

By introducing the matrix (depending on the coordinates) 

^-i = (l+i/,)/_i/,0/.^^, (15) 

these formulas can be rewritten 

n = in"^7"\ (16) 

and 

{^'^(x),^;3(y)}xO=j;0 = (7(x))„/3 5(f -y). (17) 

For the free field the matrix 7(x) reduces to 7^^. However, it appears that the Heisenberg equations, 
which are generally equivalent to the equations of motion, are not satisfied here under their most 
usual forms. To write the formulas in a condensed form, let us introduce the two operators 



V =hc 
=hc 



(1 + i/i)(-i7'5^ + mc/h) + ^/i'^S^a^ + '-(d.h^^" - a^/i)7, 
(i7'^ + mc/h){l + \h) - '-Kli^h^' - -Aduh^''' - d^hh^ 



One then finds that the Hamiltonian density 11^' — C can be written as 
while the field equation reads 



The relation 1)^ = ^1) and the formula (7"^)^ = 7'^7~^7'^ then give the conjugate equation 

Using an integration by parts, the total Hamiltonian H can be transformed into 

{d^x)'^V^ + — / ((i3x)^'(ao7"^)^'. 
The anticommutation relation ( p!7| ) then gives the commutator 



or, given the field equation. 



Similarly, one has 



then 



-ihcdoU + -hcU'ydo^ ^. 



(18) 
(19) 

(20) 
(21) 
(22) 
(23) 

(24) 
(25) 
(26) 
(27) 



The second terms in the right members of (25), (26) and (^) are unusual, since the fields ^, ^ and 
H behave as if, considered as functions of some fundamental dynamical variables, they were also 
explicitly dependent on the time. Prom the expression of H, it is in fact obvious that, among ^ and 
H one of them, at least, explicitly depends on time. But this is not obvious for 



3.1.2 Asymmetrical Lagrangian. 



The preceding difficulty takes another form if, as it is more usual |1S], the Lagrangian (|^ is replaced 
by the asymmetrical one 



C' = C + ^d,f 



The field equations are left unchanged, and the conjugate field is again: 



(28) 



(29) 



Since 11' = 11, the anticommutation relation of ^' and ^' is identical to JT^), but the Hamiltonian is 
now 



It follows that the usual Heisenberg equations are satisfied. In fact, one has 
[H', ^] = -7^^ = -ihcdo"^ , 



then 



[H',U'] 



-ihcdoU' . 



(30) 

(31) 
(32) 

(33) 



The equations satisfied by ^ and 11' are those of variables having no explicit time dependence, which 
is the rule for canonical variables. On the contrary, expression ( p9[ ) shows that depends explicitly 
on time since 7~^ does, and equation (32) is in agreement with this dependence. 

However, the Hamiltonian H' is not Hermitian. In fact, since one has (^'2?^')^^ = \I'P^', it 
follows that 



H 



't 



(d^x)^P^ = F' + i^c / (d^x)^(ao7"^)^ 



(34) 



The asymmetry of the equations for ^ and ^' is a consequence of this lack of Hermiticity. 



3.1.3 Trouble with the interaction picture. 

A common shortcoming of the two preceding methods of quantization is the absence of a coherent 
definition of the interaction picture. In fact, if this picture was defined, there would be a unitary 
operator U, such that the corresponding field variables ip and w would be given by 

^ = [/^fC/-i, w = UUU-^. (35) 

Moreover, these variables would be free fields so that one would have the relation w = ihipjQ. Such 
a relation is not compatible with ( ^5|) since one has H = ih"^j^^, 7^^ 7^ 70, and since U must be 
unitary. 

In what follows, the quantization is defined in a way such that the preceding difficulties do 
not appear. In particular, the interaction picture will be defined, allowing the construction of the 
transition probabilities and that of the S matrix. 



3.2 A coherent method of quantization 



The afore-mentioned problems will be solved by a change of variables eliminating the derivative 
coupling. The Lagrangian density @ can be written under the form 



(36) 



Let us introduce the new field Q by the formula ^ = AQ, where A is a matrix to be determined, 
which depends on the coordinates. One has 

^7-^50^) = 67° At 707-iA(5oe) + 67° At 707-1(30 A)e. (37) 

The terms containing the time derivatives of G and © in (|3^) will be those of the free Dirac 
Lagrangian if one has 'y^A'^^^^~^A = 7" or, assuming the inversibility of A, if 
By writing the matrix A under the form A = MU where M is Hermitian and positive [^] and U 
unitary, the preceding equation becomes 



7O7-1 = (M-i)2 = / + i/i - i/i°'^7°7/.- 



(38) 



If /i^jy is sufficiently small, this matrix is inversible and defines a unique matrix M positive-definite 
[34|, while U may be arbitrary. However, to ensure that the asymptotic states deduced from G or 
^ have the same physical interpretation, it is necessary that U goes to the identity when t ±00. 
More precisely, we will impose the condition U = I, by which, when h^i, = 0, the field is a free 
field as ^. 

The Lagrangian density takes the following form as a function of 0: 



C = he 



- (e7°(9o0) - (9o0)7°e) + - (erHdke) - (5fc0)r^0j + -qtq 



1^ 



(39) 



A, 



(40) 



where and F are defined by 

= j^A^° \{i + i/i)7^' - i/i'^'^' 

r = i (^j^A-^doA + T^A-^dkA^ 

- ^7° ((aoAt)At-i7° + (5,At)At-i(70rS0)) 7° " 2^(1 + lh)j'A^l'^ • (41) 
Let us note the Hermiticity relations 

The field equations, equivalent to those of ^ and ^, are now 



(42) 



i (7°(9o0) + T\dkQ)) + 2 + iidkT")) = 0, 

i f(ao0)7° + {dke)r'^) + fr - iidkF")) = o. 



while the current reads 



or, more explicitly^. 



f = c07°At7O [(1 + 1/1)7^ - Ihf^'j^] A0 , 



f = c07^0 



0^ 



(43) 
(44) 

(45) 
(46) 



^The expression of is identical to that of the free-field case. In Appendix B this property is taken as a condition 
allowing the introduction of the field Q in the framework of the linearized theory of General Relativity. 



The conjugate momentum Hq of Q has the same form as in the free-field case 



dC — 

It fohows that the anticommutation relation of G and Q is the usual one 

„o 



(47) 



(48) 



It is equivalent to the anticommutation relation (^) of ^ and The Hamiltonian density TIq now 
reads 



ihc 



he 



ere, 



which gives for the total Hamiltonian 



or, equivalently, 



He = -Y I (d^^) i ( eV'idkQ) - (afce)r^-G ) + ere 



He = -he J (d='x)e iv'^idkO) + ^ (r + ^(5^?'=)) e 



he / (cfx) 



1^ 



e. 



(49) 
(50) 

(51) 
(52) 



This operator is Hermitian, and from the preceding expressions and, from the field equations (^) 
and (H), one checks the Heisenberg equations 



[He,e] 



ihedoQ 



ihedo@ 



[i70,ne] = -ihedoUe ■ 



(53) 



The difficulties discussed in Section 3.1 have disappeared. In particular, the basic variables 0, G 
and He look like variables having no explicit dependence on time. On the contrary, the variables 
^, ^ and n, initially considered, depend explicitly on time, since one has 



(54) 



and since the matrix A generally depends on time. 

Let us remark, however, that the Hamiltonian Hq differs from the Hamiltonian H introduced at 
the beginning in terms of It is convenient to consider H from (|23|) as the integral of the density 



and Hq from (^l|) as the integral of the density 

He = hee 



r + iidkT' 



e. 



(55) 



(56) 



A rather tedious transformation of the former expression, using the expressions of 7 ^, F^, and 
the definition of 0, leads to the formula 



n^=ne + ^67° [A-HdoA) - {doA^)A^-' 



e. 



Conversely, one has 



ne = n^ + ^^7° [At-i(9oA-i) - (5oAt-i)A-i 



(57) 



(58) 



With the choice made above of the matrix U, one has A = M = Mt, with M defined by (|8|). The 
two preceding formulas can then be written 

n^ = ne + —97° [M-\doM] e, (59) 

We = + ^^7° [M-\doM-^] ^. (60) 

It is noticeable that the equality Tiqj = Tio, is valid at first order with respect to the h^^^s. In fact, 
since d^M and doM~^ are first-order quantities, this approximation is obtained by taking, for the 
first term in the commutators in ( |59| ) and (|60|), the zeroth-order approximation of M~^, that is the 
matrix unity. 



4 Interaction picture and the 5*— matrix 

With the Lagrangian (|39|) the interaction picture is easily defined, since the corresponding conjugate 
momentum He has the same form as that of the free theory. This allows us to define the evolution 
operator in that picture, then the S'-matrix. 



4.1 Evolution operator and transition amplitudes 

The field in the interaction picture will be denoted by 0. Let us recall that this operator is obtained 
from the Heisenberg operator 9 by a unitary transformation such that the field equation becomes 
the free one [^] . Accordingly, the evolution equation of the state vector reads 

ih^^m)) = Hr{t)\m), (61) 



where, in the absence of derivative coupling, which is the case for the Lagrangian (pj9|), the Hamil- 
tonian Hjit) is equal to the interaction Hamiltonian expressed in terms of 9. From the expression 



(50) of the total Hamiltonian one gets 

Hi{t) = J{d^x)niM, (62) 

with 

Wint = -'^oiT'^ - I'me) + '^{d,e){^^ - ^')e - ^(t + 2^)^. (63) 

The evolution operator in the interaction picture is, from (^), the solution of the following 



equation together with the initial condition |21] : 



ihjU{tM) = Hi{t)U{tM), Uito,to) = I. (64) 



The perturbation theory is then obtained from the integral equation, equivalent to (|64|), 

U{t, to) = 1-1 f Hj{t)U{t, to)dT . (65) 
n J to 

In what follows, we are interested in the transition amplitudes to first order with respect to the 
h^yS. These amplitudes will be obtained from the first-order approximation with respect to the 
Hamiltonian Hi of the f/ -operator, namely 



[/«(t,io) = -^ f Hi{T)dT. (66) 



to 



By introducing the Hamiltonian density (|63| ) in normal form^, and the initial and final states of the 
transition, we will consider the amplitudes 



{<^>f\u^^)it,to)m = -^l\ 



(67) 



where the pair of colons : : denotes as usual the normal product. 

The first-order expressions needed for the evaluation of (|67|) are now derived from ( ^8|) which 
gives 



7 7/., 



then, from (|^) and (|l|). 



The corresponding expression for the Hamiltonian density Timt is 



(68) 

(69) 
(70) 



n 



int 



+ 



ihc- 

1 

4 

ihc 



(dke) 



7m 



(dkO) 
7> 



+ 



he- 



me 



^ /iOO--(afc/iO^)(7S°7/.-7M7V) 



(71) 



In Section 5, the expression of the Hamiltonian Hj{t) which will be used is the expression obtained 
from the space integral of ( |7l| ) by performing the integration by parts of the term containing dkO, 
which gives 



Hi 



(t) = J {d^x)e{x)j''vG{x)e{x), 



(72) 



where the operator Vcix), acting on 0(x), is given by 

he 



Vg{x) 



-7 



^/.o° + ^a,S-7°(7V-7^7'=) 



ihc n 

H 7 

2 ' 



2/1°%° + h^'lj - /i"'7 



00„,fc 



dk ■ 



(73) 



This form of the Hamiltonian is closely related to the equation of motion (0^) of the Heisenberg 



field G. In fact, with the help of ( pOj) and (70), one checks that this equation can be written: 

[ihc-f^df, - mc^) G = 7° VgG. (74) 



In what follows, the initial and final states, which appear in (|67|), are some one-particle or 
antiparticle states. They are defined from positive or negative energy solutions Xi s-^id X/ of free 
Dirac equation 



{i-f'^d^-mc/h)xk = 0, Xk{il^%i + mc/h) = k = i,f. 



(75) 



^See To first order this prescription suppresses an infinite contribution, due to tfie energy of tfie vacuum, 

in tfie transition amplitudes for ttie antiparticles only. It can be seen, therefore, as an expression of the symmetry 
between particles and antiparticles. 



Denoting by \xk) the corresponding states, one has then to calculate the matrix element 



{xf\■.eix)J'VGix)eix)■.\x^). (76) 

As a free field, the operator 6 can be written 

e{x) = J2 / (d'p) [K{p)X^^J{x) + 4(p)xy (x)] , (77) 

r=l 

where br{p) and dr{p) are the annihilation operators for the particles or antiparticles, respectively, 
and x^p^r positive or negative energy solutions of the free Dirac equation given by jl^, ^ 



= tA^J^^^'^P)^-''^'-''''^'^'^^^ (79) 



with E{p) = c^/p^ + m^c^, p = \\p\\. In terms of these, any solution x with positive or negative 
energy can be written 

X{x) = Y,j {dWih^){^S ^X) , (80) 
the scalar product of two solutions being defined by 

(X1,X2) = j {d^x)—i{x)i\2{x) . (81) 
From (^0|) we have the following expression of the one-particle states 

IX) = E / id'p){xt\x)bl{p) \<Po) , (82) 

r 

where |0o) is the vacuum state, the correspondence x \x) preserving the scalar product. Denoting 
by 6^^^ the positive frequency part of 6, this last formula implies 



O'^+Hx) Ix) = E / id'p)x'^ri^)MP) \X) 

r 

= Y.j(dW^r\^)ixir\x) 



= x{x)\<t>o)- (83) 

The matrix element (1761) reduces to 



{Xf\eM{x)J'VG{x)e(+\x)\x^) = Xf{xh''VG{x)x^{x), (84) 
so that, from (|72|), the amplitude (|67|) reads 

{Xf\U^'Ht,to)\x^) = j\d^x)JJ{x)^'^VG{x)x^{x). (85) 
This formula will be analyzed in more detail in Section 5. 



4.2 S*— matrix 



The S'-matrix can be defined when the external field vanishes at the limit of infinite tim^. It is 
obtained from the U operator by taking the limits tQ — > — oo and t — > +00, and the corresponding 
amplitudes can then be given from (85). However, it is convenient here to return to the formula (|6^) 
in which the Hamiltonian density is given by (^)- In fact, we want to show that this expression 
can be transformed into a covariant one, the expression already given by Gupta [T^ . 

The initial and final states being the same as in the calculation leading to (p5[), formula (|6^) 
with the expression (|7l|) of Wint gives, for the first-order 5-matrix, 



{Xf\S^''^\x^) = -^ {d'x)\-Xf 



mc , 



,00 



Xi 



(86) 



This expression can be simplified by integrating by parts the term containing the derivative dkhP^- 
This leads to the two terms 



^h^' [(5fcX7)(7S°7M - 7/.7°7')x. + X7(7S°7m " l,l'l'){^kX^) 
In each of these, by introducing either of the formulas 



7S°7m - 7^7%'' 



7m7°7'] 



= 2[<5y-5y + 7V7M], 
one can insert the derivatives in the combinations ^^^d^Xi or d^Xfl^j yielding 

^h^'xI^^K - K)Xi + '^h^' [x7(5° - 7m7°)(7'5/cX.) - idkXfl')iS' " 7°7m)x. 
Adding this contribution to the remaining terms in (|86|), one gets 



(87) 



Xi 



{d'x) \ i^h^'xfXi + -^h^^XflM - ^k)x^ 



Xfl^l'^in^dkXi) - {dkXfl^h'^l^iXi 



Finally, using the Dirac equation in the last bracket, this expression reduces to 

(X/I S^^^ \xr) = -i f {d''x)'-h>^'xf (iX - Xl^) Xi- 



(89) 



(90) 



This formula agrees with the rule given by Gupta : the expression appearing under the integral 
sign is, up to the sign, obtained from the interaction Lagrangian £ — £9 in (^ by replacing ^ by 
Xi and \I' by x70. This result is identical with the one valid in the case of a nonderivative coupling, 
where the interaction picture exists directly. More generally, the validity of the rule asserted by 
Gupta, implying the use of the interaction Lagrangian defined by (^) and the covariant form of the 
propagators, can be proved at higher orders in the quantized theory defined in Section 3. 

'^Here, this condition is realized, for example, in the case of a gravitational wave. 

*The contribution of the term {h/2)Co of (W) vanishes in the first order considered here, since Xi and Xf are solutions 
of the free Dirac equation. But this property is limited to the first order. 



5 Calculation of the relativistic phase shifts in the weak— field ap- 
proximation 



In this final section we use the tools and the material derived in the previous sections to make an 
explicit calculation of the various contributions to a gravitationally induced phase shift in matter- 
wave interferometry. We restrict ourselves to one-particle or one-antiparticle states. The application 
of the formalism to many-particle states and coherent beams of massive particles of different spins 
will be developed in another publication. We assume that the incoming particles or antiparticles are 
described by the state vector |x(*o)) = \Xi) ^■t some time to before interaction (to can be conveniently 
taken to be — oo). At a later time t, this state evolves into |x(0) which interferes with a reference 
beam described by |xref) with which it is recombined in the final beam splitter. In practice, |xref) 
is produced by the other arm of the interferometer and in many cases, one will have |xref) = \Xi) ■ 
We are thus interested in the spinorial wave function for one~(anti)particle states: 

x(x) = {MO{x)\x{t)) (91) 

where 6{x), |x(0)i I'Po) are respectively the free-field operator, the one-(anti) particle and the 
vacuum state vectors in the interaction representation. It is easily shown that, when pair creations 
are neglected, this expression is equivalent to the Heisenberg field amplitude (0,in| 0(x) |<I>). 
The interference signal itself is given by the projection: 



/ 



d^xxli{x)x{x) = d^x{xrc{\0Hx)\<Po){MO{x)\x{t)) 



{Xre{\U{t,to)\x{to)) (92) 



where the space integral is over some detection volume. More generally one should consider a 
detection hypersurface a{x) and the projection: 



dcT^XrciixHxix) = / da^ (xrcfl Oix) |0o) 7m {M e{x) \x{t)) . (93) 
In this paper, we shall limit ourselves to the calculation of the amplitude: 

(Xref I X(i)) = (Xref I f/(t, to) |x(io)) • (94) 

and take the phase of this complex amplitude as the phase contribution of the perturbing interaction. 
The resulting spinor (^) will also be derived using two different methods: first, in configuration 
space and second, in the momentum representation. 

5.1 Calculation in configuration space 

The evolution equation of the state vector in the interaction picture is 

^^Jt '^^^^^ = ^'^^ ^^^^^ ' ^^^^ 

where the Hamiltonian Hj{t) is 

Hi{t) = [ Sx e\x)VG{x)e{x) , (96) 



and where the operator Vg{x), acting on the field operator 0{x), is given to first order by ( [73| 

Vg = ^mc^O/.oo + ^4S(7V-7^7') 

Z o 



ihc n 



+ ^^'Si - /i°V dk (97) 



that we shall write: 



Vg{x) = A{x) + ^-^djB^{x) + ihB^{x)dj = A{x) + \ {ihdj,B^{x)]_ 



(98) 



with: 



A{x) 
B\x) 



2^ 



Ik 



(99) 



Prom equation (91) we check that the evolution of the one- (anti) particle spinor is governed by the 
equation: 



ihdtx = -ihc^^^^djX + mc^^^x + ^GX 



(100) 



to which we may add terms corresponding to diagonal magnetic dipole and off-diagonal electric 
dipole interactions |]l^, 12 1. This equation has been used in references [|lO| and [12| to discuss all the 
terms that lead to a phase shift in an interferometer. 

To obtain the corresponding amplitude, we can start directly from the integral form of (95): 



\x{t)) = \x{to)) 



he 



d'^x'9\x')VGix')9{x') \x{t') 



(101) 



to 



So that, to first order: 

iXreilX^'Ht)) 



{Xrei\U'^'\t,to)\Xi) 

__i f d^x'xU^')VG{x')xi{x') 



to 



I 

he 

t 



d'^x { xi 



to 



d x{x 



to 



^ 

Lref 



(■ref 



A{x) + ^-^djB^x) + ihB^{x)dj 



Xi 



V/i°° + V- h] - a + ^V-h 



i ^ ^ ^ 1 
+-E • V X /i - - 
4 4 



which follows also from Gupta's form. We have used the definitions: 



1 



(102) 



a = 7*^7 



a 
a 



a 
a 



h = {h^''}, h = {h'^} . 



(103) 



The calculation of the spinor itself, by a stationary phase method in configuration space, is outlined 
in Appendix C. 

In equation ( |102| ), the first three terms lead to the familiar phase shifts of the Linet-Tourrenc 
formula [^3|, the fourth term is the spin-rotation interaction and the last two terms ensure her- 
miticity. But, because of the Dirac matrices, the interpretation of the various terms in configuration 
space is not so transparent and, in previous works non-relativistic limits have been taken either 
directly as in our own work |12] or through a Foldy-Wouthuysen transformation as in |13] in the 
case of inertial fields, to put equations ( |100| ) and (|102| ) in a form where the significance of the terms 
is more obvious. 

In this paper, we may as well take advantage of the flat Minkowski space-time and, therefore, 
we will use rather the momentum representation in the following. As we shall see, the interpretation 
of the terms is then much easier, even in their relativistic form. 



5.2 Calculation in the momentum representation 

The free-field operator 9 is written as before as: 

(+)< 



r=l 



e{x) = J2 I d^P hr{p)x'^^J{x)+dl{p)x^7j{ 



(104) 



where hr{'p) and ^^(p) are the annihilation operators for the particles or ant ip articles, respectively, 
and x^-^'* are the positive or negative energy solutions of the free Dirac equation given by (|7^ ) and 

Let us introduce the Fourier transforms B^{k,t), h^^{k): 



A{x) 
h^v{x) 



1 



(27r)3/2 
1 

(27r)3/2 

1 

(27r)3/2 

1 



d^k&{k, t)e*^-^ 
d^khf,^{k, t)e^^-^ 
d^kVG{k,p,t)e'''-^e'P-^/'' 



(27r)3/2 

VG{k,p,t) = A{k,t) - B{k,t) ■ (p+ ^hk) . 



(105) 
(106) 
(107) 
(108) 
(109) 



Let us expand x{x)j Xrcf{x) and Xii^) iii plane waves using the expansions of 9^^\x) and 9^^\x) 
For particles, the output spinor is: 

X{x) = {cPo\O^^Kx)\x{t)) 



r=l 



I d'pbr{p)x^^J{x)\xit)) 



with |lp,r) = bl{p)\cl)o) and with, using again (101), 

(Xrcflx(i)) = {Xrc{\U{t,to)\Xi) 



(110) 



(XrcflxO-^ d'2;'xLf(^')VG(x')Xi(xO 

(Xref IXi) - T- 



(111) 



n (27rn)3 



^ dt' y d^x' I d^p d^p'iXrcl I Ip,, 



S(p)'' ''''^ (2^)3/2 

7x(-')(p-')e^(p"-'-^(p')*')/^(lp.,,|Xi) 



ti^^A; VG(fe,P,t')e'''-"' 



V ^(p') 

(XrcflxO-T^ / '^i' / d'p 



(112) 



to 



d3A: 
(27r)3/2 



(X 



rcf 



p+hk,r 



I E{p + hk) V E{p) 



nWt(p + M)VG(fe,P,t')'"^'''^(p) 



(113) 



and similar expressions for antiparticles. In many cases \xi) and |xref) can be conveniently taken as plane waves, but it is usually more interesting to 
consider wave packets. Replacing |xref) by |lp',r-)) we check that the momentum representation (Ip-,. |x(t)) of x(x) satisfies: 



ihdtiMbrmxit)) 



E{p) Y E{p- hk) 



(114) 



This equation is, in momentum representation, the analogous of equation (100) in configuration space and we shall give its first-order solution later. 
It leads to a discrete set of coupled equations for a fixed or negligible recoil momentum. We illustrate below, in the case of the scattering amplitude, 
how the matrix element, which appears in the second member of (113) and (114) can be evaluated. 
To first order, the scattering amplitude (113) is: 



i X /■* f f d^k 



P+f*'-']! E{p + hk)^ E{p) 
nWt(p + nfc)VG(fc,P,t')^^'^'^(p)e<^(P^+'^^)-^(^'^]*'/'^(l-^, Ix,) . 



(115) 



The next step takes benefit from the smallness of hk/mc or of hkc/E(p) to expand the various quantities in this expression to first order in these 
parameters. The energy E{p + hk) can be expanded in a Taylor series: 



E{p + hk) = E{p) + ^^+^^^ + ... = E{p) + hk.v + M + ..., 
where 5 is the recoil shift. Prom the general transformation law of spinors in Lorentz boosts (see for example [26]): 



u{p) 



cosh 



(f) 



+ n ■ a sinh 



(f)].(0), 



where n is the unit vector along p, and tanh(^ = |p| /{'jmc), we derive the following infinitesimal transformation for spinors: 



I, n ■ a , n ■ dp 1., 

u{p + op) = < 1 H — tanh (p h - smh (p 

z p 2 



a ■ 6p 
P 



+ i sinh' 



(I) 



n X 



P 



eU(p) 



(116) 



(117) 



(118) 



we have to first order in hk/mc: 



u'{p + hk) = u'^ip){l + 



1 hk 
27 mc 



■ a 



V Kk 

X— -S}, 



2(7 + 1) mc mc 



(119) 



where 7 = Xj \J\ — /?^, = fey + 7^±, the indices || and _L designate vector parts respectively parallel and perpendicular to p. The term proportional 
to d represents a boost (velocity change) and the term proportional to S a rotation (Thomas precession). 
Rather than to calculate directly 



E{v + hk) V ^iP) 



it is simpler to calculate first the matrix element: 



E{p+ \hk) y E{p- ^hk) 
^2 



{p+ \hk)VG{k,P- khk,t)u^''\p - \hk) 



A{k,t) - B{k,t) -p 



u^'''\p- \hk) 



mc (,)^ 



E{p) 



A{k,t) - B{k,t) ■ p 



u^'''\p) + ^^u^^'^^ip) { — ■ a,A{k,t) - B{k,t) ■ p} u^'-'^p) 



2' 

1 mc? 



A^E{p) 



{r')i 



mc 



mc 



4(7 + l)E(p) 

where {A, designate (anti)commutators. The first line gives: 



^Mt(^ <! — X — • %A{k,t) -Bik,t)-p) u^^'^p) , 
mc mc 



,(-')/ 



mc (,)^ 



E{p) 



vS'>\p) 



A{k,t) - B{k,t)-p 

,2 



u^''\p) 



5rr' -^-^{kxh)- w^'''>\a± + 7CT||)u;(''') 
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(120) 



(121) 



(122) 



(123) 



where w^^^ are Pauli two-component spinors corresponding either to helicity eigenvalues or to the two values of the z-component of the spin in the 
rest frame, and where 



d= (ct_l + 7'?||) 



(124) 



is the spatial part of the Thomas-Pauh-Lubanski 4-vector operator [30, 31]. The second hne gives the term: 



^||.«t(^-^ 1^ • d,Aik,t) - ~Bik,t)-p\ uir')if) = xi.p). .«ta.(^') 



The last line gives the Thomas precession terms: 



4(7 + 1) E{p) 
ihhoo 



mc mc 



+ 



4m(7 + 1 

The last line can be rewritten to yield: 



X p) ■ w^^"aw 



2m(7 + 1) E{p) 



'-^ikxp)-w^^^''aw^^'^ + 



ih 



4:E^{p)m{'y + 1) 



{k X p) ■ h ■ p 



p ■ '^aw^ 



E{p+ \hk) V E{p- ^hk) 



^\ ^ u^'^^{p + \hk)VG{k,p- ^hk,t)u^'-'\p- \hk) 



ih 



2m(7 + 1) 

If we replace now p by p + hk/2 in order to calculate (115) this introduces the additional terms: 



{k X p) 



c^p'^hf^^.p" 
2^2 (p) 



hc^k ■ p I r c^p- h ■ P 



AE{p) 



E^{p) 



hk { -> h • pc 



If we introduce the 4-vector k'*: 



kPc = k ■ V 

K, = k 

which corresponds to the energy-momentum 4-vector exchanged during the interaction, these terms can be rewritten: 

"2 I' - ,hK^c Hk^^' 



Our final result is thus, for the scattering amplitude: 
(x„t|C""(t,to)lx(«o)) = 



——I 5^^/ 



+ 



ih 



2m(7 + 1 
ihc 



{k X p) 



2E'^{p) 



(r)]-> (r') 
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(132) 



To obtain the outgoing spinor, one can replace (xref I by (^o| G{x) in the previous expression, which gives this spinor as a sum of outgoing plane-wave 
spinors: 



to 



(27r)3/2 V E{p) 



P^h^u (p" - hK") 1 + 



2m) 



5rr' + 



ih 



2m(7 + 1) 



{k X p) 



c^p^h^^p" 
2E^p) 



ihc 
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A; X i h — h ■ 



pc 

W)J\ 



(133) 



in which an explicit phase factor is associated with each outgoing plane wave component and which can also be obtained directly from the first-order 
solution of equation (114). 

One can also perform the calculation presented in Appendix D, which gives this spinor in the form of Dirac matrices multiplying the initial spinor 
plane wave components 



2E{p) 



(p'^ + ;t/t'')Vp"U- 



hn c 



ih 



^ik.x^i[E(p+nk)-E(p)]{t'-t)/h^ 



E{p) 



(r') 



(134) 



In this formula, the integral over p can be calculated by assuming that the initial wave packet has a very narrow width in momentum space around 
a central value pq. The initial wave packet Xi{x) can then be factorized. If this approximation is not sufficient an expansion of the wave packet 
around pq can be used [22] . The k integral can be performed by turning each term linear in k into a spatial derivative and the result of Appendix C 
is recovered for the spinor in configuration space. Finally the time integral can be worked out in many cases and expresses energy conservation [22] . 



For the scattering amplitude, the comparison of equation (132) with equation (102) shows new terms directly related to the momentum exchange: 
a generalized Thomas precession and a generalized spin-gravitation interaction. To illustrate how this phase shift calculation is done from equation 
(132) we shall rewrite this equation without the terms that obviously do not contribute to the phase and use expression (116) for the energy difference, 
in which we neglect the recoil shift 6 



I / (c^'^f)«rcf(p)"i(p)/?r,refA 



+ : 



cVvv(xo+a^^0p'- 



X p 



he 
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pc \ 



2m(7 + 1; 

V X \^h{xo + vt', t') - h{xo + vt', t') ■ j 
where we also made explicit the centers of the wave packets and their polarization: 



(r)t- (r') 



(135) 



,ref 



(136) 



(the assumption that the wave packet is broad enough for a*^^{p + hk) ~ a*^f{p) has been made). 

If we assume, for simplicity, that the reference wave packet is identical to the unperturbed wave packet: a^^f (p)/3r,ref = 0(i{p)Pr,i and that they are 
very narrow in momentum space around a central value p, the phase simplifies to: 



Sip 



t ( 
dt' ' 



to 



\2E{p)' 



p^h^^{xQ + vt\t')p^ + 



7 



m(7 + 1) 



c^p^'Vh^y{xQ + vt',t')p^ 
2E^{p) 



X p 



■ s 



V X h{xo + vt', t') - h{xo + vt', t' 



pc 
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(137) 



□ 



where s is the mean spin vector^ 

1*= ^p*^-(3r^^ihw^'-^^aw^'''^/2-f (139) 

r,r' 

In fact the phase calculation is usually more involved since the previous formula applies only 
to the case of straight unperturbed trajectories. In practice however, one cannot always ignore the 
fact that, when calculating the phase to first-order for a given term of the Hamiltonian, the motion 
of the particles is affected by other terms. One example is the calculation of the gravitational shift 
within the atom beam splitters, in which one cannot ignore the important effects of the diffracting 



electromagnetic field on the trajectories of the particles |3g, |3^, 40 1. Gravitational phase shifts 
have to be calculated along these trajectories. Another example is the gravity field itself, which, 
on earth, gives parabolic trajectories for atoms. The phase shift for the other terms in equation 



(137) has to be calculated along these parabolas. A convenient way to achieve these calculations is 



to replace xq + vt' and v in equation (|137D by the classical trajectory {x{t'),v{t')} obtained in the 



ABCD formalism developed in references 37|. 



Expression (137) displays all the terms which may lead to a gravitational phase shift in a matter- 



wave interferometer. They are summarized in Table 1 where one finds successively: 

• the terms involving /iqo lead to the gravitational shift (/iqo = —2 g ■ r/c^), to shifts involving 
higher derivatives of the gravitational potential and to the analog of the Thomas precession 
(spin-orbit coupling corrected by the Thomas factor). 

• the terms which involve h = {h^''}, give the Sagnac effect in a rotating frame {h = Q x r/c), 
the spin-rotation coupling and a relativistic correction (analogous to the Thomas term for 
/loo)- They describe also the Lense-Thirring effects coming from inertial frame-dragging by a 
massive rotating body, which is a source for h. 

• the other terms, which involve the tensor h = {h^-'} describe genuine General Relativity 
effects such as the effect of gravitational waves and de Sitter geodetic precession (which also 
includes the Thomas term for /loo)- 

Our expressions are valid for spins and 1/2 and may be conjectured to be valid for arbitrary 
spin if a/2 is replaced by the corresponding spin operator S. 

The reader will find calculations of the phases corresponding to these various terms in references 
|, 0, 0, m, 0> il> H, El- In these calculations, one should never forget that the external field 



/i^jy acts not only on the atoms but also on other components of the experiments, such as mirrors 
and laser beams and that, depending on the chosen gauge, additional contributions may enter in 
the final expression of the phase which should, of course, be gauge independent. 

In present experiments on the earth gravity measurements, the relative sensitivity is Sg/g ~ 
3.10~^ after 60 seconds and the absolute accuracy 5.10~^||6|, For rotations, the best sensitivity 
achieved up to now is 6.10~^'^rad.s^^Hz~^/^|^^ but these numbers are expected to improve rapidly 
in the near future, especially in space experiments, in which general relativistic effects should be- 
come detectable. An accurate measurement of the effect of gravitation and inertia on antimatter 



also appears as a possibility already discussed in reference [47| with a transmission-grating interfer- 
ometer, although we believe that an antiatom interferometer using laser beams for the antihydrogen 
beam splitters (Ramsey-Borde interferometers) would be better suited for such an experiment. The 
formalism introduced in this paper to deal with antiatoms should be useful to discuss such exper- 
iments, especially when coherent beams of antihydrogen will be produced either by Bose-Einstein 
condensation and/or by stimulated bosonic amplification. 

^More generally, if we use ai{p,r) instead of ai{p)(3r,i, the mean spin vector should be written: 

fcl ^ j{d\)e\x)mx) \xr) = J id'p)aUp,r)a,ip,r')hw^^^^aw^^'^ /2j (138) 



Corresponding energy term V 


h 

"■/if 


IVjiTTiP of flip pfrpcf 




iNewtonian potential: n,oo = ^ty /c = —2,g • x/c 


Gravitational red shift 


or acceleration field /iQo = 2a - x/c? 


Acceleration shift 


Gravity gradient g{z) ■ x = — (g — g'z/2) z 




or curvature RoiojX^x^ 




Fermi gauge: h^Q = h^{t — z/c).{x'^ — y'^)/'2, + hx{t — z/c)xy 


gravitational waves 




1 'Y r ^ 1 
/ioo = 2U/(? gives V = — VU x p -s 

mc^ 7 + 1 L J 


Thomas precession 


—cp - h 


r? ni" p "t" in o" ivi^iwf^' n — (/ V i* 1 0'\'\t(^*^ T/ — — \/ • // 

J- WJ Li Cli 1)111 ci 11 CXlllC' • 1 V u ti /\ Ju 1 \^ til V Co \ J u ' ■ 


kJClgllCi*^ C-llC-V^lj 


hf\- P'lvpn nv f np TiPTmp — T'nirrTnP' mpf.rip 

\j% O C-Ll '-'j^ ljll\Z> J_JdlO\j J- 1111 1 111^ lllO Ul 1^ 


T.pTiGp — T'TiiTrinP' (orKifpil j 

■ in 1 1 r^\. _I_ ±±±± X lllt^ 1 KJj. Ul LCtl 1 


-(c/2) 


V X h ■ s 


h = ^x x/c gives V = —Vt ■ s 


Spin-rotation interaction 




hoi given by the Lense-Thirring metric 


Lense-Thirring (spin) 


7 

m(7 + 1) 


V 


(^cp- h/E^ X p ■ s 




~ Thomas for rotation 


cV- h .p/2E 


Schwarzschild metric in isotropic coordinates: 

hm = hii = /i22 = /i33 = 21/ /c^ gives 
V = p^U/E in addition to EU/c^ from /iqo 




Einstein gauge: hn = — /122 = h+{t — z/c), 
hi2 = h2i = hxit- z/c) 


Effect of gravitational 
waves 


(c/2) 


V X ^ 


T • pc/E^ 


■ s 


Schwarzschild metric: U = —GM/r 

hoO = ^11 = /i22 = /i33 = 211 /c^ 

gives V = — K— VC/ X p ■ s in addition 
mc^ 7 L J 

1 T r ^ Jl - 

to K = — K V(7 X p ■ s irom /loo 

mc^ 7 + 1 L J 


de Sitter 
or geodetic precession 




Einstein gauge: 
hij 


Interaction of the spin 
with gravitational waves 


7 

2m(7 + 1) 




^c^p- ^ -p/E^ 


j X p • s 




~ Thomas for gravitation 



Table 1: Classification of the various energy terms entering the expression of the phase shift. The factor 7 is the time dilation factor and should not 
be confused with the PPN parameter ^yppN which can also be introduced (for de Sitter precession it gives the familiar factor {^ppN + 1/2))- 



5.3 Analogy with the electromagnetic interaction 

The formulas that we have derived for the transition amplitude and for the ongoing spinor strongly 
suggest analogies with the electromagnetic field case. To emphasize these analogies, let us introduce 
the following pseudo-potential 4-vector 



(143) 



then the Linet-Tourrenc term, which appears also in the generalized Thomas precession is simply 

.2 1 . 



(144) 



where ii^ is the 4-velocity p^/m. 
The corresponding field 



{E/c,B) 



(145) 



appears in the outgoing spinor (|134D through 



and in the generalized spin-gravitation interaction term 



ihc 
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k X \ h — h ■ — 7— 
E{p)) 



2^E{p) 
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(146) 



This new correspondence between the gravitational interaction and the electromagnetic interaction 
generalizes the so-called gravitoelectric and gravitomagnetic interactions introduced by de Witt [45| 
and Papini |4(:]. 



Appendix A: Dirac equation in curved space— time 

For a given space-time manifold, together with its metric tensor g^^, the Lagrangian density of the 
Dirac field reads [^J] 



he , — — 
2 V y 



il^T^ a — rnc/ti 



he , — — 



-iV — mc/h 
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with 



V — P'-t^ 



V, 



(148) 



In these formulas the eg's are four- vector fields constituting an orthonormal tetrad (tetrad indices 
have a hat), that is satisfying the relations 



(149) 



More rigorously, one should introduce [p3| 



which stems directly from a compact form of the interaction Hamiltonian 

with p° = —a-'pj + ■y^'mc and pj = ihd 



(140) 

(141) 
(142) 



with {r]ai3) = diag{+l, —1, —1, —1). A change of tetrad is possible according to the formula 

e^'^ = A^„e^^ (150) 

where A is a matrix of the Lorentz group defined at each point of the spacetime manifold. In such a 
change the Lagrangian density (147) is left invariant if the Dirac field is correspondingly transformed 
according to the law 



^r' = 5(A)-1^', 



(151) 



the matrix S{A) being the usual transformation matrix under a Lorentz transformation of a Dirac 
spinor [0]. For an infinitesimal transformation 



Aq,/3 = rja^ + eal3, 



one has 



5(A)-i = / + ^e,^a"^. 



The field equations read 

ij°'^a — mc/h 



^' = 0, ^ 



-iV a'y°' — mc/h 



0. 



(152) 



(153) 



(154) 



With the invariance of C under the phase transformations ^ — > ^'e ^' — > ^'e*°, is associated the 
current density 



(155) 



where J'^ is a four- vector invariant under a change of tetrad, and one has the conservation relation 
in either one of the two equivalent forms 



(156) 



Appendix B: Weak-Field approximation 

It is now assumed that space-time admits a coordinate system (xf^) in which the metric tensor takes 
the form 



livi 



(157) 



According to that hypothesis, the h^^s will be considered as first-order quantities, and the subse- 
quent calculations will be valid at this order. To determine the corresponding form of the Lagrangian, 
it suffices to construct the tetrads associated with ( |157D . By putting es^ = rfa^ + /q^, where fa^ is 
of the first order, one obtains from (149) and ( |157| ), 



/a/3 + //3o — hafj ■ 

The general solution of these equations is of the form 



(158) 



(159) 



in which the e^^'s are first-order quantities only restricted by the antisymmetry condition e^p 
—epa- One then has 



(160) 



Introducing the free Dirac Lagrangian (|2|) and the associated stress-energy tensor (P), the Lagrangian 
density ( |147| ) calculated at the first order and corresponding to ( p. 601 ) reads^ 

£ = £o - k^^'T^u + '^e^^'^i^^t, - ^5,7^)^ " ^d^e.p-^i^a'''' + a^^^)^. (161) 
2 o 

The corresponding field equations are the following 

+ '-{d^h - a,/i%,)7^ + ie'^'^^X - \d^e,p^^'a''P]^ , (162) 

= ^{-iJ^^'^^,-mc/h){l + lh) + ^"^u-ff,h^"' 

-^{d^h - d,h\)j^' - i'a ,7/.^^" - Id^^upa'^Pj''] ■ (163) 

The terms depending on h^u are identical to those appearing in (^) and (^). The same is true for 
the current density, which is now 

f = c^[^M + i/i^A* _ i/i/^-^, _ e'^-'-f^]^. (164) 



Considered independently of the context, the equations (|162| ) and ( |163D are invariant under the 
transformations of the Poincare group provided that the h^^^s and the e^^^'s are transformed like the 
components of a second-rank tensor, and \I' by the corresponding transformation law of a spinor. 
However, the weak-field character of h^i, is not conserved by any finite Lorentz transformation 
except by the rotations. Moreover, the /i^,y's being the basic quantities, the condition that the e^,y's 
are of first order is naturally interpreted by assuming that these latter are some linear functions of 
the former, that is one has 



a^,p„hP\ (165) 



Such a relation is not compatible with the general tensor transformation law, on account of the 
symmetric or antisymmetric character of h^u or e^y, but can be made compatible with the rotations 
by a suitable choice of the coefficients Of^^pa . 

The rotational invariance of (165) is equivalent to the relations 

^/.V'p'a' ^ r/ R^'^' r/ R^^' aP'^P" , (166) 

where R can be any rotation matrix. The corresponding solutions are given byP| 

^OiOj ^ i^^ij ^ ^Oijk ^ Q ^ ^jkoi ^ i^^ijk ^ ^ijki ^ Q ^ (-^g7) 

where A and B are some arbitrary parameters. If, in addition, the invariance under parity is 
postulated, one has to take B = 0, giving 

eoi = -eio = -Ahoi , eij = , (168) 

then, instead of ( |16[1| ), 

% = 1 + 1^00 , % = (l - A)hoi , ejg = (i + A)hoi , err. = rjij + ^hij . (169) 
In particular the choice A = corresponds to the tetrad 

^^The indices of ft^^ and e^^ are raised with the help of 77'"'. 

^ As usual, the latin indices can take the values 1,2 or 3, and e'-' is the completely antisymmetric symbol. 



which we cah the standard tetrad associated with h^^. The corresponding Lagrangian is identical 
to the Lagrangian (|l|), therefore the equations (||) and @ are recovered from (162) and (163) by 
letting e^jy = 0. We will continue to designate the corresponding field by ^. 

If we write explicitly the spinorial connection in this weak-field approximation 



(171) 



the equation for ^ can also be written in a simple form, analogous to the electromagnetic case, and 
which will find an interpretation in terms of gauge fields in the flat space-time approach of the main 
text 



Z7 



0. 



(172) 



Another choice of tetrad can be made in relation with the expression of the vector current. From 
155|) written in the weak-field case, one has in general 



and then 



f/c = (1 + i/i)^7^^ - i/i^''^'7^^ - e^'''^!-i^^, 



f/c = (1 + \h)^-fH - ah^" + e°'^)^7i.^- 



This quantity will be proportional to the usual density ^'7"^' if one has e*^* = - 
obtained with the choice A = ^. Denoting by the corresponding field, one has[^ 



(173) 



(174) 



-i/i'^*, which is 



3 



c(l + i/i%)^'7°^' • 



It is then possible to introduce the new field Q by |12] 

G = (l + i/i^)5^'~(l + i/i%)^', 

such that 

f = ce7°e, 



(175) 



(176) 



(177) 



as in the free-field case. The change of field \I' ^ ^I'' corresponds to a change of tetrad as defined 
above, so that, from ( |160| ) and (|170| ), one has the formulas (151), (152) and (153) written with an 
infinitesimal parameter e equal to — e that is such that Sij = 0, eoi = — £^0 = ^/loi- This gives 



then, from (176), 



^' = (/-i/ioi7V)^, 



e = (/ + i/i-iV7V)^- 



(178) 



(179) 



The field is identical to the field introduced in Section 3.2 for the purpose of defining the quanti- 
zation and the interaction picture. 



^^Let us note that, in the expression of the conserved charge corresponding to (175), the volume element {l + ^h\)cfix 
is, in the first-order approximation, the 3-volume associated with the spatial metric defined from the coordinate system 



Appendix C: A stationary phase calculation 

The outgoing spinor can be calculated directly in configuration space with the help of a stationary 
phase formula. In the particle case, from the definition 

X{x) = {MO^^\x)\x{t)) (180) 

and the integral equation 

\x{t)) = \x{to)) - ^ l\d^x')e^+^\x')VGix')e^+Hx')\xit')) (isi) 

one can derive, to first order, the expression 



X{x) = Xr{x) - J- I ((i^x')(-/.o|e(+)(x)^(+)(x')|0o)7°VG(x')(<Ao|e(+^(x')|x^), (182) 



'to 

which, with the help of standard formulas, can be transformed into 

ft 



{x)=x^{x)-^ f <fx'S^+\x-x')^'^VG{x')x^{x) (183) 



= ^*(^) -yJ^I "^'P^^ri^) I {d''x')x'^^J{x'h^VG{x')x^{x'). (184) 
By introducing (|85|), one finds 

Xix) = X^ix) - ^^Yl [ dW^H^) iMx^^r^\U^'Ht,to)\x^). (185) 

r 

In this last expression the matrix element can be submitted to the same transformations as those 
leading from ( |8^ ) to (|0[), yielding 



X{x) = Xi{x)-iY, dW^^\x) d^x'^h^'^ix')x^+J{x'){^^X-dllf?jXiix') 

= X^{x)-\ f d^x'h'^''{x')S^+\x-x')(jX-Zl,)xr{x'). (186) 



By introducing the expression 



one finally gets for a plane wave with the momentum pQ 



po=+E{p)/c 



(187) 



X{X) = X^ix)-'-^ fdt' /dV/l^-(x')x 



to 



(2^ / |^(^' + ^'o).(7^Pp + "^c)e-^(---')/S^x.(x'). (188) 



The elementary stationary phase formula^^ appUed successively to the integral over p, then to the integral over x', yields the Linet-Tourrenc 
result [23], namely : 



to 



Xi{x) , Vq 



Pqc 

Em 



This expression is the beginning of an asymptotic expansion in powers of h of the form 

„2 rt 



^~2hJt dt'F{x-vo{t-t'),t') 



Xiix), 



(189) 



(190) 



in which F is defined by an expansion in integer powers of h whose first term corresponds to the Linet-Tourrenc formula. This expansion can be 
obtained from a generalization of the stationary phase formula given by Hormander [25] . The application of this general formula to the two integrals 
appearing in (188) yields, after some complicated calculations, the following expression of the outgoing spinor (to be derived in the next Appendix 
by a simpler method) 



E{p) 2E{p) 



E{p) 



E{p) 



5'^ 



c c 



didj {p^Puh^ 



(191) 



v=po 



Appendix D: Derivation of the wave function using the momentum representation 

The method used for the amplitude in the main text can be used to derive the wave function. The spinorial wave function for one-(anti)particle 
states is 



For particles: 



(</.o|0W(x)|x(t)) = {Ue^+\x)\xi 



X{x) = (0o| 0{x) \x{t)) ■ 



(192) 



h (2vr/.)3/2 



to 



mc 



{2Tr f E{p + hk) 



—u^''\p + hk)u^''^^{p + hk) 



VG{k,p,t')e^'-'eii'^'^^+'^^^^^^^ 

V E{p) 



(193) 



*See, for instance [24]. 



n 



□ 



and a similar formula for antiparticles. To first order we get 



Vo(*,P,tOe*'"eil'^»+'*'>-'^®l(''-*Y|^«<'''®e*»--'^»){lp>.|Xi) (194) 

The next idea is to express the propagator which appears in this equation in order to write the outgoing spinor in the form of Dirac matrices 
multiplying the initial spinor plane wave components. 
From 

J2 u^'^ m^'^ iP) = 2^ [7% + mc] (195) 

r 

one obtains 



E{p + hk) 

from which, after some algebra, one finds the spinor 



^[E(j!+nk)-E{p)\{t'-t)/n ^ ^ik-v{t'-t) 
Finally one can check that: 



2E{p) I \E(p) 



e 



(196) 



where the time-dependent exponential can be expanded to any desired order for recoil shift corrections: 

[1 + i5{t' - t)] . (198) 



(Xrefl U^'\t,to) \xito)) = / d'x xU^)x^'\x) ■ (199) 



We obtain indeed 



with 







47 



from which equation (132) is recovered. 



\ m). 
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